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The measurement of azimuthal distributions in nucleus-nucleus collisions relies 
upon the assumption that azimuthal correlations between particles result solely 
from their correlation with the reaction plane (i.e. flow). We show that at SPS 
energies, the ansatz is no longer valid, and two-particle correlations due to mo- 
mentum conservation, final state interactions or resonance decays become of the 
same order as those arising from flow. This leads us to introduce new methods to 
analyse collective flow, based on a cumulant expansion, which enable to extract 
smaller values than those accessible to the standard analysis. 



In a collision between two heavy ions, the azimuthal distribution of outgo- 
ing particles with respect to the reaction plane, i.e. collective flow, is expected 
to reveal new insights on the central region of the collision: thermal equilib- 
rium, equation of state, time evolution.EI It is therefore important to have 
reliable flow values, and to be able to measure small values. 

We first discuss in Sec. |l| the method usually used to analyze flow data, 
and we point out that it is based on an assumption which is not valid at ul- 
trarelativistic energies. In the next two sections, we introduce a new method 
which allows flow measurements at such energies: the main ideas are pre- 
sented in Sec. H, while the practical implementation in terms of the event flow 
vector, including acceptance corrections, is given in Sec. |^. Finally, we briefly 
summarize our results in Sec. IJ. 

1 Flow and two-particle correlations 

The standard method for analyzing ±he experimental data relies on a study of 
two-particle azimuthal correlations.0 The one-particle azimuthal distribution 
is first expanded into a Fourier series, whose coefficients Vn characterize the 



flowfl 



(e"^), (1) 
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where (j) denotes the azimuth with respect to the reaction plane, while the 
average is performed over a large number of events. 

Since the actual orientation of the reaction plane is not known, only 
differences between particle azimuths can be measured. Thus, flow is ex- 
tracted from the measured two-particle correlations (e™^'^!^'^^)), under the 
assumption that they are only due to flow, or at least that the contribution 
^gm(0i-02)^^ of direct, nonflow correlations is negligible. In other words, it is 
assumed that in the decomposition 



the second term in the right-hand side is much smaller than the first. j-. 

However, we have shown that this ansatz is not valid at SPS energiesH 
nonflow two-particle correlations, which can stem from various sources, are of 
the same magnitude as two-particle correlations due to flow. Thus, the NA49 
values for pion and proton flow are significantly modified when two-particle 
azimuthal correlations arising from momentum conservation, HBT quantum 
effects, and resonance decays are subtracted from the measured correlations. 

Two-particle correlations may not be easy to take into account in the stan- 
dard analysis, even if the sources are well identified (as for example in the case 
of resonance decays). Moreover, the possibility that there are still unknown 
sources cannot be discarded. Thus, the safest approach consists in using a 
method which does not require a precise knowledge of these correlations. 

2 Measuring flow using multiparticle azimuthal correlations 

We recently proposed a new method for the analysis of flow, which does not 
depend on two-particle correlations since these are eliminatad by means of a 
cumulant expansion of multiparticle azimuthal correlations .0 In this section, 
we show how this method allows measurements of integrated flow (i.e. aver- 
aged over a phase space region, see as well as differential flow for given 
transverse momentum and rapidity ( |2.2| ) . 

2. 1 Integrated flow 

The main idea underlying our method is the following. In a collision where 
N particles are emitted, direct /c-particle correlations are typically of order 
so that they become smaller when k increases. In particular, their 
magnitude decreases faster than the flow contribution to the measured k- 
particle correlations: one can expect that the comparison between both terms 
allows measurements of smaller and smaller flow values. 
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Let us consider for instance the measured four-particle azimuthal cor- 
relations. These correlations can be expanded as the sum of products of 
direct one-, two- and three-particle correlations, where "one-particle correla- 
tion" means correlation with the reaction plane, that is flow. The order of 
magnitude of each term in the expansion can easily be estimated, keeping in 
mind that a-flow term is of order w„ while direct fc-particle correlations are 
0(1/A^'^^^).EI After calculation, the dominant terms in the expansion are: 

[For simplicity, we have assumed that w„ is not much smaller than V2n, and 
neglected a term 0(-y|„/A^^).] The first and the fourth terms in r.-h. s. are 
of magnitude and respectively: a direct comparison between them 

would allow flow measurements down to values Vn ^ i /N^^^, smaller than the 
^ 1 /iV^/^ limit of the standard analysis which can be deduced from Eq. (H). 
In order to compare these two terms, we have to get rid of the second and third 
terms, which are in fact equal, and correspond to two-particle correlations. 

That can be done, taking the cumulant of the four-particle azimuthal 
correlation. The purpose of this cumulant is precisely to remove lower order 
(i.e. involving k' particles with fc' < fc = 4) correlations. It is defined as the 
difference between measured correlations (e'"("^i+<^2-<^3-<^*)) - 2(e™('^i-'*='))2. 
Using Eqs. (|) and (|), one finds that the cumulant is: 

^^g««(0i+02-03-04)\^\^ = _^ ^g»n(0l+02-03-04)^ = ^ O ^ (4) 

so that it gives the flow, provided this latter is larger than Measuring 
this cumulant is indeed sensitive to smaller flow values than the standard 
method. For instance, at SPS energies where about N = 2500 particles are 
emitted, using the cumulant allows measurements of w„ ^ 0.3%, while 
the standard method is limited to Vn ^ 2%; thesa^alues are to be compared 
with the published values ui ~ t;2 — 3% for pions.El 

2.2 Differential flow 

Let us assume that the average value Vn — (e""^) of flow is known for a given 
particle species, which we shall call "pions" to fix ideas. For any type of 
particles, let us say "protons", although it can be anything — even "pions" — , 
it is possible to perform detailed measurements of the flow v'^{pT, y) = (e™"^) 
using multiparticle correlations in order to go beyond the limitations of the 
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standard two-particle methods To siraplify the expressions, we shall assume 
n = m; however, the ideas remain valid is m is a multiple of n. 

Let us consider for example the four-particle azimuthal correlation be- 
tween a proton and three pions (e™*^'^^'^^^'^^^'^^^). As above, this measured 
correlation can be expanded in terms of lower order direct correlations, yield- 
ing a term v'^^v^, a term of order corresponding to direct four-particle 
correlations, and other terms. These latter, and in particular two-particle 
correlations, can be removed taking the cumulant [by analogy with Eq. (^] 

Inspection of this equation reveals that the measurement of the cumulant in 
the left-hand side gives access to as soon as it is larger than l/{Nvn)^, 
which is also the accuracy on w^. 

3 Practical implementation using the event flow vector 

3. 1 Integrated flow 

An easy way to implement the method outlined in section ^ consists in using 
the event flow vector, which is defined as 



i=i 

where the sum runs over particles from the same event, and M should be 
chosen as large as possible. Averaging Eq. (|^) over many events, one finds 
{Qn) = VMvn'- a nonvanishing (Q„) signals collective flow. 

From definition (^, it is obvious that the powers of \Qn\^ involve mul- 
tiparticlc azimuthal correlations, and only differences between angles (which 
require no knowledge of the reaction plane orientation). That explains why 
the method can naturally be expressed in terms of the fiow vector: the average 
value {Qn) corresponds to w„, i.e. what we want to extract from the data, while 
the measured quantity (IQnl'') corresponds to measured fc-particle azimuthal 
correlations. 

Let us turn once more to four-particle correlations. By analogy with 
what was found in section |^, the expansion of the moment (|(3„|'*) involves 
various terms: (Qn)^, a term (|(5n|'*)c corresponding to direct four-particle 
correlations, and other terms which we wish to get rid of. These latter can 
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indeed be removed, using the cumulant of the Qn distribution, which at order 
fc = 4 is defined as 

{{\Q.X)) ^ {\Qn\') - 2{\Q^\Y = -{Q„)^ + 

Since {Qn) ~ VMvn, this equation is strictly equivalent to Eq. (^), provided 
M and N are of the same magnitude. Therefore, the measurement of the 
cumulant (Q) provides a way to extract smaller Vn values than what can be 
obtained within the standard two-particle analysis. 



3.2 Acceptance corrections 

An import£LDt feature of the method is the easy implementation of acceptance 
corrections The only modification regards the definition of the cumulant: 
in Eq. the cumulant involves only two moments of the Qn distribution 
because the other terms vanish for a perfectly isotropic detector. If the detec- 
tor is not isotropic, these terms should be reintroduced in the definition, and 
the flow vector Qn measured in the laboratory frame is not equivalent to the 
flow vector Qn with respect to the reaction plane. Taking into account these 
changes, the cumulant is (dropping the n for brevity): 

mn^W) - 2 (Q) {QQ*') - 2 (Q*) {Q*Q') - 2 im' - {Q') {Q*') 
+8 (Q) (g*) (m + 4 (g)' (g*2) + 4 (q*)' (q^) - e (Q)' (Q*)' 
= -(g>* + o(i/Af). (8) 

This expression generalizes Eq. (|^) to the case of real, nonperfect detectors, 
and should be used to extract flow from experimental data. 



3. 3 Differential flow 



Let us readopt the same notations as in Sec. 2.1. Representing the pions' 



azimuths by the flow vector Qn [Eq. (o)], the four-particle (1 proton, 3 pions) 



correlations considered in Sec. 2A appears in the moment {\Qn\'^Qn^^"^)- The 
idea is again to expand this moment, and to take the cumulant in order to 
remove unwanted terms. More precisely, the corresponding cumulant for the 
(l+3)-particle azimuthal correlation is: 

((|g„pQ;e'"'^)) ^ {\Qn\^Q*y-^) -2{Q*y-^) (|g„p) 

= -(g„)X-Ho(^^). (9) 

Once the pion flow (Qn) has been measured, this cumulant gives v'^ with 
accuracy 0(1/Af(g„)'^) = 0(l/(Mw„)'^), as did the cumulant (||). 
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4 Summary 



We have shown that the standard flow analysis, which rehes on two-particle 
azimuthal correlations, is in trouble when w„ becomes of order l/N^/"^. In 
opposition, the new method we advocate, based on a cumulant expansion of 
four-particle correlations, allows integrated flow measurements down to values 
Vn > l/N^/'^. It is even possible to extract smaller values, using higher order 
cumulants, as for example the cumulant of the six-particle azimuthal corre- 
lation, or equivalently the cumulant ((IQnl^)): we have derived a generating 
equation which gives the cumulant of the Qn disti'ibution at a given order, 
and relates it to the corresponding power of {Qn)^ 

This new method also allows the measurement of differential flow, from 
the correlation between a particle azimuth and the event flow vector Qn- 



Once again, the resultS|X)f Sec. ^ can be generalized to arbitrary order, using 
a generating equation.^. Finally, corrections for detector inefficiencies are 
naturally implemented in the method, using a redefinition of the cumulants. 

We wish to emphasize that all the cumulants, either of multiparticle az- 
imuthal correlations or of the Qn vector, can be obtained from measured 
quantities, that is, the moments of the Qn distribution: our method relies 
only on measurable quantities, and therefore can easily be implemented. 
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